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Abstract. We interpret the Central Limit Theorem as a fixed point theorem 
for a certain operator, and consider the prob lem of li nearizing this operator. In 
classical as well as in free probability theory | VDN92 |, we consider two methods 
giving such a linearization, and interpret the result as a weak form of the CLT. 
In the classical case the analysis involves dilation operators; in the free case 
more general composition operators appear. 



1. Introduction 

The ubiquity of the normal distribution as indicated by the Central Limit The- 
orem (CLT) is a somewhat mysterious result. One of the possible explanations 
for it is an interpretation of the CLT as a fixed-point theorem; for a plethora of 
approaches see ||Tro59| , |Gol76| , |HW84| , Par86| , |Swa9 1| , |Sin92|| . The starting point 
for this analysis is the following weak form of the CLT for independent identically 
distributed random variables. Let the operator T be defined on probability mea- 
sures on R by T\i = (u * u) o Sj_ . Here * is the convolution, and S r (for "scaling" ) 

is the dilation operator, d(/j, o S r )(x) = d[i(r~ 1 x). We call this operator T the 
central limit operator. The theorem follows from the CLT, and is well known. 

Theorem. The fixed points of the operator T are the scaled normal distributions 
X ° S t - If fi is a probability measure with zero mean and unit variance, then the 
iterations T n fi weakly converge to \- 

In this approach, the starting point for the study of the CLT is the investigation 
of the operator T . This operator is clearly non-linear, and as the first approxima- 
tion we consider the properties of the linearization of this operator [|Sin92|1 . These 
properties are of course well-known (although we have not found adequate refer- 
ences; but see ||Sin92|| and also |[5in76|| ). However, there is now a different version of 
probability theory, with its own CLT, which has not been investigated to the same 
degree. This is the free probability theory of Voiculescu (for an introduction, see 
which, in particular, turns out to describe the behavior of certain 



e.g 



VDN92 



large random matrices. In this theory the notion of (commutative) independence 
is replaced by the notion of free independence (for the operator-theoretic definition 
and the motivation behind it see e.g. ||VDJN92|| ). Now the classical convolution 
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can be denned in terms of independence as follows: fi * v is the distribution of the 
sum of two independent random variables with distributions fi and v (and it is a 
theorem that the distribution of the sum of independent random variables depends 
only on the distributions of the summands). Correspondingly, in free probability 
theory one defines the (additive) free convolution of measures fi EB v as the distri- 
bution of the sum of two freely independent random variables with distributions 
[L and v [ Voi85 , Voi86| , |BV93 , |VDN92|| (and the above comment applies). Thus 



the free central limit operator is T(fi) = (/i EB /i) o Sj_. One of the main technical 

differences between classical and free theories of probability is that the operator 
of convolution with a given measure is linear, while the operator of taking a free 
convolution with a given measure is highly non-linear. However, the above opera- 
tor T is non-linear even in the classical case, and thus one can expect similarities 
between linearizations of the classical and free versions of this operator. 

We have two somewhat different approaches at our disposal. The original one, 
initiated and largely developed by Voiculescu |[Voi86| , |BV93| , |VDN92|| (see also 



|Maa92| ), is to define a certain operation, called the R-transform, on the space 



of analytic functions, which linearizes the additive free convolution. Thus this 
operation is an analogue of the logarithm of the Fourier transform in the classical 
case. Another approach, due to Speicher ||Spe90|| and developed, among others, 



by Speicher and Nica [|Spe94| , |Nic95|| , is to use a certain analogue of the classical 



combinatorial moment- cumulant formula. This approach is somewhat less general, 
but the parallel with the classical situation is more explicit. 

In what follows we want to indicate the parallels between the classical and the 
free case. Therefore, whenever appropriate, we will use the same notation for both 
cases. The situations where (important) differences between the two theories arise 
will also be indicated. 

Acknowledgments. We would like to thank Prof. D.-V. Voiculescu for suggesting 
the problem as well as many helpful discussions. We would also like to thank Prof. 
N.G. Makarov for some suggestions. 

2. The Combinatorial Approach 

2.1. Notation. Let /x be a probability measure. We denote by T both the central 
limit operator T(/z) = (/i * //) o Sj_ and the free central limit operator T(/i) = 

(/i Sfi)oSj_. It will be clear from the context which one is meant. We denote by T 

the manifestations of T on auxiliary spaces: in Section 0, the spaces of sequences; 
in Section |^, the spaces of continuous functions. Precise definitions will be given 
at appropriate times. Also, we denote by x the appropriate normal distributions: 
standard Gaussian in the classical context and the standard Wigner semicircle law 
Voi85|| in the free context. 



2.2. Background. For a measure fi, its n-th moment is m£ = J x n dfi(x). In this 
section we consider only probability measures whose moments of all orders are 
finite. In fact, throughout most of the section we disregard the non- uniqueness 
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and identify the measure with its collection of moments. Thus let A4 be the space 
of all one-sided real-valued sequences. We will call the elements of Ai the moment 
sequences and denote them by m := {mj}^, even though only some of them 
are moment sequences of measures. On M. we set up the topology of entrywise 
convergence; this is the weak*-topology on the space M. as the dual of the space 
of the "eventually 0" sequences, and it turns M. into a topological vector space. 
Note that if a sequence of elements of M. do in fact correspond to measures, and 
if its limit corresponds to a unique measure, then one has weak convergence of the 



corresponding measures [pur 91 



For every such moment sequence m there is also the corresponding free cumulant 
sequence c := {ci}°l 



Spe94| , [Nic95|| determined by 



(1) 



(free) 



m k 



={B ly ..,B n } 



Here V nc (k) is the set of noncrossing partitions of the set {1, . . . , k} [ |Kre72| , |Spe90| , 
Spe94| , [Nic95|| , which can be described as follows: these are partitions of the ver- 
tices of an /c-gon such that the vertices in each class can be connected by lines 
inside the fc-gon so that the lines for different classes do not cross. Also, B^-s 
denote the classes of the partition ir, and \Bi\ denotes the number of elements of 
B,. 

The classical cumulant sequence can also be described by a similar formula 
Shi96| , N ic95|| , namely 



(2) 



(classical) 



m k 



c\ Bj \ 



--{B u ...,B n } 



where V(k) is the collection of all partitions of {1, . . . , k}. 

Let us denote the transformation from the moment sequence to the cumulant 
sequence determined by the formula (P (resp., (Q)) by 1Z : {m i \°^ l — > {c i }^ 1 . We 
call 1Z~ X the cumulant-moment transform. Note that TZ is given implicitly; there 
are also explicit formulas [ Nic95 |. Note also that fc-th moment depends only on 
the cumulants of orders less than or equal to k, and vice versa. Therefore both TZ 
and 7Z~ l are continuous bijections M. — >• Ai; however, we will think of the domain 
of TZ as moment sequences and of its range as cumulant sequences. 

The point of the transformation from the moment to the cumulant sequence is 
that for sequences corresponding to probability measures, the appropriate action 
of the operator T on the cumulant side is linear. Indeed, 



and 



r k 4 



where c 71 are the cumulants of the measure rj. That is, 



(3) 



2 1 - 



2^ 
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Thus define, on the space of cumulant sequences, the operator by (7^(c))fc = 
2 1 ~ fc / 2 c/ c , and on A4 the operator T = 'Rr l oT n o 1Z. Clearly, since T 7 ^ is linear, 
in order to linearize the operator T, we are interested in the linearization of the 
cumulant-moment transform Tlr x . 

In the sequel, by a linearization of a map A at a point x we mean its Gateaux 
derivative: (D x A)(y) = lim e ^ A ^ x+£y ^~ A ^ when the limit exists in the appropriate 
topology. Note also that from (|3|), a fixed point of T for which the cumulant 
sequence is defined must have all the cumulants other than the second one equal 
to 0. In the classical case, this describes the Gaussian distributions; in the free 
case, this describes the free normal distributions, which are the dilations of the 



Wigner semicircle law ||Spe90 



2.3. Proposition. The linearization of the cumulant-moment transform at (the 
cumulant series corresponding to) the normal distribution \ [respectively, standard 
Gaussian in the classical case and standard Wigner semicircular distribution in 



the free case ||VDN92|1 ) is given by a (formal) infinite lower-triangular matrix A = 
(aij)i°j=i, where 

a. In the classical case, a n+ 2k,n is (n — 1)! times the number of partitions of 
(n + 2k) elements into classes exactly one of which contains n elements and 
the remaining k classes are pairs. 

b. In the free case, a n+2 k,n is the number of noncrossing partitions of (n + 2k) ele- 
ments into classes exactly one of which contains n elements and the remaining 
k classes are pairs. 



In both cases aij = if j > i or (i—j) is odd. For explicit values, see Theorem pT^ 



Proof. The value of the n-th cumulant is a polynomial function of the first n 
moments only, and vice versa. Thus in the topology of entrywise convergence, the 
differentials of both TZ and 1Z~ 1 exist. 

Given two moment sequences m°, m d , we define the sequence {f(m°,m d ) n }'^ =l 
recursively by 

(4) <= e Efn^iW^vi 

7r={-Bi,...,B n } 

where {c°} = lZ(m°) are the free cumulants. Then 



e n fei + ^( m °' m< vi) = < + e < + 

7T = {B 1 ,...,B n } 

Note that if m° = m M , m d = m^, then the last expression above is just m M+£i/ + 
o(e). So the sequence {/(m°, m d ) n }^L 1 is the derivative of the moment-cumulant 
transform TZ at m° in the direction m d . For m° = m x , the free standard normal 



LINEARIZATION OF THE CLT OPERATOR 



5 



(semicircular) distribution, cf = 5i2, and so (HI) becomes 



k 



(5) m d k = J2ak,nf(m a 



n=l 



where f(m d ) := f(m x ,m d ) and a n +2fc,n is the number of noncrossing partitions 
of (n + 2k) elements into classes exactly one of which has n elements and the 
remaining k classes are pairs. Thus m = Af, where A is the lower-triangular 
matrix (A) it j = a^j. 

In the classical case, we start with 

(6) m d k = (i^|-l)!^(n c l^l)/ (m0 ' m<i)|Bil 

TreP(fc) i=l / 

w={B lr .. ,B n } 

and by the same sort of reasoning see that the derivative of 1Z at the standard 
Gaussian is given by the lower-triangular matrix A with a n+ 2k,n — (n — l)!x (the 
number of partitions of (n + 2k) elements into classes exactly one of which contains 
n elements and the remaining k classes are pairs). □ 

As stated above, the operator T n is linear. It is easy to see that its spectrum is 
discrete. Its eigenvectors are the cumulant sequences £j = {Sij}^, for j = 1, 2, . . . , 
with corresponding eigenvalues 2 1-J / 2 . Therefore for the linearization of operator 
T, the eigenvalues are the same, and the eigenvectors are the moment sequences 
Cj = {aij}^, where Ojj-s are defined in the above theorem. In fact, these are true 
moment sequences, and so give the eigenfunctions for the central limit operator 
T. 

2.4. Theorem. On the space of measures with all moments finite, the lineariza- 
tion of the operator T has eigenvalues 2 1 ~ n ^ 2 ,n = 1,2,.... The corresponding 
eigenfunctions are absolutely continuous with respect to the Lebesgue measure, 
with densities: 

a. In the classical case, ^■e~ x2/ ' 2 = e~ x2 ^ 2 H n (x) , multiples of the Hermite poly- 
nomials |Sin92jl . 



b. In the free case, l[^2,2](t)-^=^T n (t/2), multiples of the Chebyshev polynomials 
of the first kind. 

Proof. In the classical case, a n +2fc,n = (n — 1)! x the number of partitions of 
(n + 2k) objects into one class of n elements and k classes of 2 elements. It is easy 
to see that a n+ 2k, n = an d ak, n = for k < n or (k — n) odd. Therefore 

for fixed n the Fourier transform (defined by YlJLo jr m i(^)' ? = / e ttx dfi(x)) of the 
n-th eigenfunction of T is 

oo 

E^Mr-^(^) n+a = -(*trexp(-t 2 /2) 
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and the sum converges absolutely. Thus the eigenf unctions are the multiples 
of Hermite polynomials -^e~ x ' 2 dx = e~ x l 2 H n (x)dx (note that these are not 
exactly what one usually means by the Hermite functions). 

In the free case, a n +2fc,n = the number of noncrossing partitions of (n + 2k) 
objects into one class of n elements and k classes of 2 elements. It has been 
calculated by Kreweras [|Kre72|l to be a n+ 2k,n = ( n+ k ) ( one uses an inductive 
argument based on the following fact: a partition 7r with a class of n elements is 
noncrossing iff each of the n intervals in the complement of this class is a union of 
complete classes of it, and ir restricted to each of these intervals is noncrossing). 
The Cauchy transform (defined by EJLo m i*~ 0+1) = / ¥?) ( i Akll65 L |VBNjg| , 
see also the next section) of the n-th eigenfunction is Ylk { n+ k k ) z~( n+2k+1 \ For 
z G C + , the series converges absolutely for \z\ > 2. Its integral is 



n + 2k — 1\ ^_ (n+2 fc) _ _ 1 (n + 2k\ ^_( n+2 k) 

n + k\ k ) ^ n + 2k V k 

k x ' k x 

f2k\ „-f2fc4-n Tk„„ 17 fJ\2 



In particular, for n — 1 we have F\(z) = — ^2kk+i( 2 k) z ^ 2k+1 " > ■ Thus F\( 
—F\(z)z — 1. Therefore F\ is related to the generating function for the Catalan 
numbers [|Rio68|| , and is in fact ~ z +f^ . Similarly 

Lemma. Forn > 1 the integral of the Cauchy transform of the n-th eigenfunction 

« -i (^) , -T. k ^Tk{ n+ k k >- {n+2k) = -\ (^P 3 ) 

Proof of the Lemma. The series converges absolutely for \z\ > 2. The proof is 

/ _ / 2 _/l \ n / _ / 2_A \ m ( _ / 2_/l \ ("+ m ) 

by induction, using the identity ( — ^ — - J ( ^— ^ — - J = ( — ^ — - J . By 
equating coefficients, we have to prove the combinatorial identity 

En fn + 2k\ m fm + 2l\ n + m fn + m + 2t 
{ _ t n + 2k \ k J m + 2l \ I J ~ n + m + 2t \ t 

k,l>0 

for m, n — 1, 2, . . . and t — 0,1, ... . But this is a particular case of the generalized 
Vandermonde (also known as Rothe) identity ( [[Rio68| , Sec. 4.5], see also |pK66| ). 

The above moments determine a unique distribution ||Dur91|| , and we can see 
directly that the n-th eigenfunction is related to the Chebyshev polynomials of 
the first kind, namely it is a scalar multiple of l[-2,2] (/)^7j=7f T n {t/2) dt, where 
T n (t) = cos _1 (n cost). □ 



2.5. Remark. Besides being eigenfunctions of the operator DT on a topological 
vector space, the above functions in fact form orthogonal bases in (smaller) Hilbert 
spaces. The Hermite functions -^e~ x I 2 = e~ x l 2 H n {x) for n = 0,1,... form 

an orthogonal basis in L 2 (e x2 ^ 2 dx), while the Chebyshev functions of the first 
kind l[_ 2) 2](£)-/2=p7n(t/2) for n = 0, 1, . . . form an orthogonal basis in the space 
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L 2 (l[_ 2! 2](t)\/4 — t 2 dt). Thus, in the classical case the operator DT is a compact 
self-adjoint operator on L 2 (e x2 ^ 2 dx), while in the free case DT is a compact self- 
adjoint operator on L 2 (l[_ 2 ,2](£)"\/4 — t 2 dt). 

We have the following interpretation of the CLT. 

2.6. Corollary. We say that a measure /i is in L 2 ((p) if /x <C (pdx and G 
L 2 (ip). Note that for all the measures in L 2 {e xl I 2 dx) or L 2 (l[_ 2)2 ](t)"\/4 — t 2 dt), 
the moments of all orders are finite. 

a. On the space of probability distributions in L 2 {e x2 l 2 dx), with mean and vari- 
ance 1, the normal distribution \ as a fixed point of the central limit operator 
T is strictly spectrally stable, that is, the differential of the operator at this 
point has the spectrum inside the unit disc. 

b. On the space of probability distributions in L 2 {1 [_ 2j 2](t)v / 4 — t 2 dt), with mean 
and variance 1, the free normal (semicircular) distribution \ a s a fixed point 
of the free central limit operator T is strictly spectrally stable. 



Proof. Classical case: ||3in92|| The Hermite functions J^e x I 2 = e x ^ 2 H n (x) for 
n — 0, 1, . . . form an orthogonal basis in L 2 (e x ' 2//2 dx). The conditions on the mo- 
ments of the distribution in the hypothesis mean that in the linear approximation 
we consider only the perturbations / with 0th, 1st, 2nd moments equal to 0. This 
means precisely that / is orthogonal to e~ x I 2 and (the densities of) the first two 
eigenf unctions of T . Since the eigenvalues of T are 2 1_n//2 , all eigenvalues for 
n > 2 are less than 1. 

Free case: the Chebyshev functions l[- 2;2 ](t) ^*_ ti T n (t/2) for n = 0, 1, . . . form 

an orthogonal basis in the space L 2 (l[_ 2)2 ](t)-\/4 — t 2 dt). Again the eigenvalues 
are less than 1 for all but the first 3 of these. But here, the hypothesis correspond 
to the orthogonality to the first 3 Chebyshev functions only if all the distributions 
considered are supported in the same interval [—2,2]. Note, however, that due to 
the results in [ BV95|| , this restriction is weaker than it appears. □ 

Thus we would expect that on this subspace the fixed point is attracting, just 
as the CLT states. 

2.7. Remark. The moments a n+ 2k,k can be calculated in a way similar to the above 
for the setting of the R q transforms [|Nic95|| , related to g-independence. However, 



there does not seem to be in that case a nice recurrence formula, and so the 
corresponding eigenf unctions are not calculated directly. Moreover, note that 
the orthogonality conditions in the above Corollary are different in the classical 
and the free cases: in the free case the inner product is given by the free normal 
(semicircular) distribution, while in the classical case it is the inverse of the normal 
distribution. Also, we are not aware of any standard interpolation between the 
Hermite polynomials and the Chebyshev polynomials of the first kind. Thus one 
would not necessarily expect to have a similar construction for the interpolations 
between free and classical cases, e.g. related to g-independence. 
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3. The Analytic Approach 

In this section we consider the problem of linearizing the operator T by analytic 
means. First we briefly go over the classical situation. 

3.1. Classical Picture. Let a G (0, 2], /3 = 2~^ a , T af i = (p*fi)oSp. For (p a = a- 
strictly stable distribution ||Shi96| , |Dur91|| (the skewness coefficient does not appear 
explicitly in the sequel and so is not included in the notation), T a (<p a ) = ip a . Then 
the differential of T a at cp a is 

DT a u = 2{y * ip a ) o Sp 
Taking the Fourier transforms, 
(7) DT^{t) = 2u(j3t)<p a (/3t) 

Also by stability 0^,((3t) = a (t). Therefore for i>(t) = h(t)0 a (t), the right-hand- 
side expression in (|7|) is 2h((3t)0 a (t). For h, on the space of continuous functions 
the eigenfunctions are h(t) = t a , a G C, Re a > or a = 0, with eigenvalues 
2 • (3 a = 2 • 2~ a / a , corresponding to v a {t) = t a a (t)- Here we use the principal 
branch of the logarithm. Now let (3 = l/y/2, i.e. cp a = tp2 — X- m this case, 
among all the eigenfunctions we can distinguish the integer values of a as follows: 
among the functions t a , the smooth ones are precisely those for a G N. Thus 
among all u a , the ones whose densities decay faster than any polynomial are just 
the u n -s, n G N [5E96 , |Dur91|| . These measures are manifestly in L 2 (e x I 2 dx). 



They are v n = J-^e x I 2 dx, with eigenvalues 2 1 n//2 , and we obtain the result of 
the previous section. 

3.2. Free Picture. In the free probability picture, the main device is the R- 
transform, introduced by Voiculescu ||Voi86| , |B V93| , |VDM92|| . Given a measure /i, 



for z G C\supp(/i) one defines the Cauchy transform (sometimes called Stieltjes 
or Borel transform) of \i by G^(z) = J ^zj-. For positive \x the Cauchy transform 
is an analytic map C + <-> C~; it has the property G^(z) = G^z). The measure 
can be reconstructed from its Cauchy transform by taking the boundary values 
ImG^i + Oi) ||Akh65| , Ch.3, Addenda and Problems], [ |Hor90i 3.1]. 



On a nontangential neighborhood of (Stolz angle) in C~, we can define 
K^{w) = G~ 1 (w), and the /^-transform R^{w) = K^{w) — ^. is an ana- 
lytic map — > C~ on a nontangential neighborhood of 0. The main prop- 
erty of the i?-transform is that it also linearizes the additive free convolution: 
R[Mv{w) = Rfj,(w) + R u (w). In fact, if all the moments of a measure \i are finite, 
then Rp,{w) = cf-u/ -1 , where Cj are the free cumulants. Also R^ s r ( w ) — 

rR^rw). Thus the action of the operator T a on the i?-transform side is just 
RT a ^{w) = 2/3Rn(/3w), and in particular it is linear. Here we again define the 
operator 

T a (fj) = (fi ffl fj) o S,3 
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for a G (0,2], and (p a = free a-strictly stable distribution ||BV93| , [Pat95| , |BPB96| . 
By the above observations, the linearization of T a is again given by the linearization 
of the i?-transform. 

The quasilinear differential equation governing the behavior of free convolution 
semigroups has first appeared in [ VoISB , Theorem 4.3]. Here we need a variant of 
that theorem. The proof is quite similar to | |Voi86| ] . 

3.3. Theorem. Let u be a freely infinitely divisible probability measure ||Voi86| , 
|V13N92|| . Let ip be a function analytic in C + . For any point z G C + ; for small 
enough t the function G(z,t), which is the functional inverse of K v (z) +tip(z), is 
defined at z. Consequently for all z G C + ; 

8C 

(8) G' u (zMG u (z)) + —(z,0) = 

Proof. By | |B V 93| , Proposition 5.12] for v freely infinitely divisible, G v maps C + 
conformally onto a domain. Let Q be a bounded domain in C + whose closure is 
contained in C + . Choose t so that for z G £1, t \ip'(G u (z))\ < 1. Then the function 
z +tift(G v (z)) is univalent on Q and invertible on its image. Denote this image by 
Q' and this inverse function by f t . Since, for z£fl, 

{K v + tip)(G u (z)) = z + tip(G v (z)) 

we also have 

{K v + ty){G v (f t (z)) = z 

for z G Q'. Consequently we can define G(z, t) = G u (f t (z)) for z G Q'. 
As one possible construction, define 

= { z\ — n < Re z < n, — < Im z < n 
I 1 n 

Let t n satisfy the condition above and also t n \ip{G u {z))\ < 1/n for z G VL n . Then 
G(z,t) is defined on 

iY = < z I — n -\ — <Re^<n , — < Im z < n 

1 I ' n n n n 

These domains exhaust C + . On their common domains, since their inverse is 
analytic in t, G(z,t) is different iable in t; the derivative exists as a limit in the 
topology of uniform convergence on compact sets in the upper half plane. 

By definition G[K V + tip(z),t) = z for z G G U (Q). Differentiating with respect 
to t at t = 0, we get 

F)G f)C 

—(K u (z) 1 0Mz) + —(K„(z),0) = 

for z G G U (C + ). Substituting G v (z) for z, we get the required equation, for 
zeC+. ' □ 
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3.4. Remark. In the above theorem no consideration is given to the positivity or 
even existence of the distribution corresponding to the Cauchy transform G(z,t). 
While we do not know of a satisfactory description of these, there are various 
conditions. 

3.4.1. Necessary Conditions. For K v + if) to correspond to a positive measure, it 
is necessary that (1) the nontangential limit of zif)(z) as z — > be 0, and (2) there 
exist a Stolz angle at 0, T C C+ s.t. (K v + if>)(T) C C" iPTOf . 



3.4.2. Sufficient Conditions. 

a. The following are the known cases where K v +ip corresponds to a positive mea- 
sure: (1) if) is an i?-transform of a positive (hence necessarily freely infinitely 



divisible) measure [Voi86|. (2) v is the free normal (semicircular) distribu- 
tion; if) is analytic in a neighborhood of the unit disc, sufficiently small, and 
if)(z) =ip{z) pV95[ . 



b. For a measure /i and a distribution u, the i?-transform of ((1 — e)n + ev) is 

R {1 - £) , +eu = R^w) - eK'^w) ■ (G v - G^K^w)) + o{e) 

Denote if) = K'Jw) ■ (G v — G ^(K^w)). Then if v is positive, the deformation 
in the direction of if) is tangent to a curve of positive measures. Note also that 
we have the inverse formula, 

G v (w) = G,(w) + V(G*(«0) • G ',M 

3.5. Discussion. The theorem can be interpreted as follows. Let G be a univalent 
function on C + . Define K = G^ 1 , 71(G) (z) = K(z) — 1/z, and let if) be analytic 
on G(C + ). Then the derivative of the map 7Z~ X (which is nothing other than 
the differential of the operation of functional inversion) at a point 7Z(G) in the 
direction if> is —G'if>(G). This linear map is invertible; the inverse linear map 
(which is the differential of TZ) is if) i— > —K'if)(K) (where if> is now analytic in C + ). 
Finally, let T^(if))(z) = 2f3ip(f3z) and denote G a := G lfia . Then the derivative of 
T a = 7Z- 1 oT^oK at G a is 

D G T a (i))(z) = 2G' a (z)PK' a (f3G a (z))i)(K a (f3G a (z))) 

= 2if)(u a (z))u' a (z) 



Here 



u a (z) = K Va ((3G Va (z)) = -u Vaa s , Va 

-I X -k r 



is a particular instance of the transition probability function of ||Voi93| , |Bia98|1 . 

The eigenfunctions of on the space of all analytic functions in the upper half 
plane are of the form te l ^z a , with eigenvalues 2 l ~ <yCl+l ^ a . Restricting to various 
spaces selects particular values of <f>, a. Thus the eigenfunctions for the differential 
of T a (resp., T a ) are the (boundary values of) the functions e % ^G' a G a a . 
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3.6. Example. For the free normal (semicircular) case a = 2, v = x the Cauchy 
transforms of the eigenfunctions of the operator DT (which are the eigenfunctions 
of the operator DT) are given by 



X X 



V7 2 



By taking the boundary values — ^ Im G(t + Oi) (see [|Akh65| , |Hor90f ) , the eigen- 
functions themselves are 



+ 



e w sin ( y log 



t- 



+ 6 + X7T 1 1 



L (-oo -2] 



(t) dt 



t- v? 



e y7T sin ( y log 



+ 



-[2,oo) 



(t) dt 



+ 



: exp(ycos 1 (t/2)) cos (x cos 1 {t/2) — <pj 1[_ 2,21 (^) dt 

\ \ I- 

It is easy to see that the Criterion |3.4.1| requires that (for some <fi) we have x = 
Re a > 1 or (x = Re a > —l,y = Ima = 0); the corresponding point spectrum is 
the union of the unit disc and the interval [1, 2]. On the other hand, for — 1 < x < 
1, y — the functions are the i?-transforms of freely stable distributions, while for 
x G N, y = the functions are entire. Thus by Criterion |3.4.2|a| the corresponding 
eigenfunctions are in the tangent space to the space of positive measures. 

All the moments of a measure are finite iff its Cauchy transform is analytic at 
infinity. If a function G is defined by the above expression on C + and satisfies 
G(z) = G(z), it is analytic at infinity iff a G N and = 0. Notice that among 
the above measures, these are precisely the compactly supported ones. Explicitly 

1 



their Cauchy transforms are 
are 

1 



The eigenfunctions themselves 



: cos in cos 



vT - 7I ^\x/2))l h2 , 2] (x)dx 

That is, we recover the Chebyshev functions of the first kind. 



3.7. Remark. Here we can see another difference from the classical case. As noted 
above, by a result of Bercovici and Voiculescu [ BV95|| the deformations in the 
directions z n ,n G N actually produce positive measures (for small enough time). 
This is in contrast with a classical theorem of Marcinkiewicz, which states that 
for P a polynomial, e p is never a characteristic function (i.e. a Fourier transform 



of a positive measure) if the degree of P is greater than 2 (see e.g. |[Ram67| , Thm 
3.13]). 



3.8. Example. For the 1-stable symmetric distribution, which is the Cauchy 
distribution (fx, the eigenfunction Cauchy transforms are , and the eigen- 
functions are D a ^-^ dx. Notice that these are exactly the same as in the classical 
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case. Therefore not only are the free 1-stable distributions the same as classical 
[B V93|1 , but their small neighborhoods look the same as well. 



ones 



3.9. Remark. For a general freely stable distribution, using for example the for- 
mula DgTZ -1 ^) = —jtt(G), one can obtain parametric expression a la Biane 
BPB96| , Appendix] for the densities of the corresponding boundary values. In 



particular, one has such expressions for the densities of the eigenfunctions of the 
free stable central limit operators. It is not clear whether they are of use. 

3.10. Composition operators. Finally, we have a brief discussion of the con- 
nections with the theory of composition operators (see e.g. | |Val31| , |C1V195| ). The 



action of the operator T a on the primitive (in C + ) of a Cauchy transform is, up to 
a multiplicative constant 2 and up to an additive constant, just the composition 
with the function u a . 

One of the main theorems about composition operators is that such an opera- 
tor is necessarily conjugate to an operator of composition with a linear function 
3M95| , Thm. 2.53]. In our case, in the terminology of | |CM95| , Section 2.4] the 



operator T a has a natural halfplane-dilation model provided by conjugating with 
the linearization of the i?-transform, i.e. G a o uj a = (3G a . A fundamental set 
CJM95| , Defn. 2.54] for yt> is an open, connected, simply connected domain A such 



that C A and the iterates of any compact set end up in it after a finite 

number of steps. It is not hard to see (e.g. |[BPB96| , Appendix]) that C + serves as 



a fundamental set for both G a and u a while G(C + ) is a fundamental set for K a 

A number of results on the spectra of composition operators on various classical 
spaces are known. In particular, on a Hardy space H°° the spectral radius of 
a composition operator is equal to 1 | CM95 , 3.1]. In our case the composition 



operator is defined on tangent spaces to a certain cone in H°°. 
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